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The exact quantum decay of a one-dimensional Bose gas with inverse-square interactions is pre¬ 
sented. The system is equivalent to a gas of particles obeying generalized exclusion statistics. We 
consider the expansion dynamics of a cloud initially confined in a harmonic trap that is suddenly 
switched off. The decay is characterized by analyzing the fidelity between the initial and the time- 
evolving states, also known as the survival probability. It exhibits early on a quadratic dependence 
on time that turns into a power-law decay, during the course of the evolution. It is shown that 
the particle number and the strength of interactions determine the power-law exponent in the lat¬ 
ter regime, as recently conjectured. The nonexponential character of the decay is linked to the 
many-particle reconstruction of the initial state from the decaying products. 
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Understanding the decay dynamics of unstable isolated 
systems is of relevance to a wide variety of fields ranging 
from quantum science [1, 2] to statistical mechanics [ 

5 ] and cosmology [ ]. While the exponential decay law 
is ubiquitous in Nature [7], quantum mechanics dictates 
its breakdown when the time of evolution is large [ ] or 
small [ ]. The existence of these deviations follows from 
the linearity of unitary quantum dynamics. Consider the 
preparation of a unstable quantum state. During its de¬ 
cay, the time-evolving state can be decomposed as a co¬ 
herent superposition of the initial state and a set of decay 
products (any state orthogonal to the initial state). Devi¬ 
ations from exponential decay result from the possibility 
for the decay products to reconstruct the initial state. 

To be precise, let I’I'o) = |^(0)) be an unstable quan¬ 
tum state prepared at t = 0, that evolves into T(t)) when 
the dynamics is generated by a Hamiltonian H(t). The 
probability to find the time-evolving state l’I'(f)) in its 
initial state I’I'o) is referred to as the survival probability 

m, 

5(t):=|A(t)| 2 = |(vI/o|vk(t))| 2 , (1) 

where A(t) = (’l'ol'I'(t)) is the survival amplitude. Equiv¬ 
alently, the survival probability is given by the expec¬ 
tation value of the projector on the intial state V = 

l^oX^ol, 

S(i) = tr[|*(i))<*(i)|P] (2) 

and is identical to the fidelity between the ini¬ 
tial state and the time-evolving state, S(t) = 
F[|4' 0 )(^ , oU^ , (i))(^ , (OI] I ]• The survival probability is 
as well closely related, but different, from the Losclimidt 
echo [ 3 ]. However the distinction is often omitted in re¬ 
cent literature. 

To appreciate that the decay dynamics of S(t) is gen¬ 
erally non exponential, it is convenient to use the Ersak 
equation, that relates the values of survival amplitude at 
different times during the course of evolution [ 9 ], 

( 3 ) 


Above, the memory term is given by 

= (^o\U(t, t)QU (t, 0)I’I'o), (4) 

where U(t,t') = Texp(— i J* H{s)ds/h) is the time evo¬ 
lution operator from time t! to t, and Q = 1 — V is 
the complement of the projector on the initial state V. 
Hence, M(t,r) represents the probability amplitude for 
the initial state to evolve into decay products at time r, 
and subsequently reconstruct the initial state at time t. 
For an exponential decay to hold, the memory term is 
to vanish. However, the memory term plays a dominant 
role at short and long times of evolution. The short-time 
decay is known to be governed by the energy fluctua¬ 
tions of the unstable state. Experimentally, it was first 
demonstrated in [11] and its existence sets the ground for 
the quantum Zeno effect [12]. It is a consequence of uni¬ 
tary time-evolution, provided that the first and second 
moments of the Hamiltonian exist, see [13, 1 ] for excep¬ 
tional cases. That deviations from exponential decay are 
as well to be expected at long times was pointed out by 
Khalfin in 1957 , for systems whose energy spectrum is 
bounded from below [< ]. Measurements consistent with 
these deviations were reported in [15]. 

The decay dynamics of many-particle quantum sys¬ 
tems has recently received a great deal of attention [16- 
26]. These studies show the need to characterize quan¬ 
tum decay at a truly many-particle level, beyond its de¬ 
scription in terms of one-body observables such as, e.g., 
the integrated density profile [27-30]. Achieving this is a 
challenging goal, due to the limitation of reliable numer¬ 
ical techniques. Even at the single-particle or mean-field 
level, propagation methods based on space discretization 
in a finite spatial domain can introduce artifacts due to 
the enhanced reflection from the boundaries of the nu¬ 
merical box, unavoidable for long-time expansions [31]. 
An attempt to palliate this effect with complex absorbing 
potentials [32] explicitly suppresses state reconstruction, 
and delays the onset of power-law behavior in an un¬ 
physical way [33]. By contrast, these issues are absent in 
studies of fidelity decay in spin systems [34]. As an out¬ 
come, analytical results in quantum decay, often based on 


A(t) = Ait — t)A{t) + M{t , r). 
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time-dependent scattering theory, are highly desirable. 
Recent theoretical progress has been mainly restricted to 
two particle systems [16-22] and quasi-free few-particle 
quantum fluids [23, 24, 26]. Experimentally, the role 
of Pauli exclusion principle has been demonstrated us¬ 
ing analogue simulation in photonic lattices [35], while 
interaction-induced particle correlations have been mea¬ 
sured in optical lattices [36]. 

In this article, we present the exact quantum decay dy¬ 
namics of an interacting many-body system that is equiv¬ 
alent to a gas of particles obeying generalized exclusion 
statistics. We rigorously show that the survival probabil¬ 
ity decays as a power law at long times with an exponent 
that depends on the strength of the interactions and the 
particle number. In the non-interacting limit, this result 
proves the scaling conjectured based on the study of few 
particles systems [16, 17, 22, 23, 25]. The non exponential 
character of the evolution is linked to the multi-particle 
reconstruction of the initial state. 


I. MODEL 


The Hamiltonian model we consider is that of N bosons 
effectively confined in one-dimension in a harmonic trap 
and interacting with each other through an inverse- 
square pairwise potential. This is the so-called Calogero- 
Sutherland (CS) model [37, 38] 


N 

m=t[ 


i—l 


h 2 d 2 1 o 2' 

2^a? + 2 m " w 


i<j 


h 2 A(A- 1) 
m\qi - q 3 \ 2 ' 


equation 

b + K(t)b = &“ 3 , (6) 

where K(t) = [u>(t)/u o] 2 , and the boundary conditions 
6(0) = 1 and 6(0) = 0 follow from the stationarity of the 
initial state. Note that SU (1,1) and the scaling dynamics 
are robust against the breakdown of integrability [46, 47]. 

The ground-state of the CS model has a Bijl-Jastrow 
form [38] 

i N 2 

*=1 j>i 

Here, the normalization constant Cm,a reads 


N—1 


C N ,A = 2-^ 1+A ( N - 1 >](2 7 r)t Yl 


3= 0 


r(i + ( J + i)A) 
r(i + A) ’ 


it is related to the normalization constant of the prob¬ 
ability distribution function for the Gaussian (/3 = 2A)- 
ensembles in random matrix theory and derived using 
Mehta’s integral [38, 48, 49]. Using the dynamics (5), it 
is found that 


( - / N,A 


5n,a(I) ^N[l+A(N- 


1 )] 


/R N 


11 ** 

i—l 


o- q -i-( 1+A- 


'III* - q 0 

j>i 


2A 


Upon explicit computation, the following closed expres¬ 
sion is obtained 


Its energy spectrum (purely phononic) and the complete 
set of eigenstates are known [37-39]. The CS Hamilto¬ 
nian includes several relevant limiting cases. It reduces 
to non-interacting bosons for A = 0, while the Tonks- 
Girardeau gas [40-43] describing hard-core bosons is re¬ 
covered for A = 1. For A ^ {0,1}, it describes interacting 
bosons which are equivalent to Haldane anyons [44]. We 
shall assume that the frequency of the trap for t < 0 is 
given by ui{t < 0) = ujq. From now on, we use dimension¬ 
less variables q —> q/qo, t —> t/to with q 0 = ( h/mujo ) 1 / 2 
and to = Wo" 1 - 

To describe the decay of the survival probability, we 
first note that the CS model belongs to a broad class 
of systems for which the exact time-dependent coherent 
states can be found [45, 46]. A stationary state 4/ of the 
system (5) at t = 0 with energy E, follows a self-similar 
evolution dictated by the SU(1, 1) dynamical symmetry 
group, 


^(qi, ■ ■ • ,qn,t) = exp 
0 2 


-qr' 0 

2=1 




( 5 ) 


where r(t) = f* dt'/b 2 (t'). Here, the scaling factor 
b = b(t) > 0 is the solution of the Ermakov differential 


S N ,A(i) = [a(f)6(t)] -N[1+A(N-1)I , (8) 

where the function a(t) is given by 



As a result of the boundary conditions, 6(0) = 1 and 
6(0) = 0, a(t) reduces to unity at t = 0. For N = 1, 
one recovers the survival probability of a single particle 
in a time-dependent harmonic trap Si.\(t) = a(t)b(t). 
It follows that the survival probability of a N particle 
CS system is identical to that of N non-interacting parti¬ 
cles obeying Generalized Exclusion Statistics (GES) with 
exclusion parameter g = A. GES was introduced by Hal¬ 
dane for systems with a finite Hilbert space [44], and 
extended by Wu to unbounded Hamiltonians [50]. It ac¬ 
counts for the number of available states excluded by a 
particle in the presence of others, and it smoothly extrap¬ 
olates between bosonic and fermionic exclusion statistics, 
and beyond. Note that the many-body-wave function 
is always symmetric under the permutation of particles, 
this is, the exchange statistics is always bosonic for arbi¬ 
trary A. The exclusion parameter is defined as the ratio 
g = — Ad/AN of the change in the available states Ad 
as the particle number is varied by AN. For the CS 
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model the exclusion parameter is precisely given by A 
[ ]. Particles with fractional excitation A 7 ^ {0,1} are 
Haldane anyons. From ( 8 ), the survival probability for 
N non-interacting and hard-core bosons is recovered for 
A = 0,1 respectively. This leads to the duality relation 

Sn,a ( t ) = [<Sn,o(^)] 1 A [Sn,i(£)] A , (10) 

that represents a signature on the quantum decay dy¬ 
namics imprinted by the transmutation of statistics ob¬ 
served in a CS system as a function of the GES param¬ 
eter A. Indeed, Eq. (10) resembles the known relation 
for the equilibrium partition functions obeyed by Haldane 
anyons [51]. More generally, 

<Sn,7+a(£) = ['SN, 7 (i)] 1 A [SN,7+i(t)] A . (11) 

Simly put, these mathematical identities emphasize the 
fact that the CS model smoothly extrapolates between 
non-interacting and hard-core bosons (or more generally, 
between different types of Haldane anyons). In addition, 
this duality is not only apparent in equilibrium proper¬ 
ties as discussed in [ ], but can be clearly manifested in 

nonequilibrium observables as well. 


II. SUDDEN EXPANSION 

Suddenly switching off the trap (e.g., K(t) = 0(— £)), 
leads to an expansion dynamics with the scaling factor 
given by b(t) = y/1 +1 2 > 0. The survival probabil¬ 
ity decays monotonically as a function of time and there 
is a smooth transition between the short and long time 
asymptotics, as shown in Fig. 1 for different values of A. 
At short-times, 5n,a(£) is a quadratic function of time, 

S*,x(t) = 1 - In 1 + A(N - 1 )}t 2 + 0(f 4 ), (12) 

O 

where all odd moments vanish identically. Given the gen¬ 
eral short-time asymptotics of the survival probability 

S n ,a(*) = l~AH 2 t 2 + 0(t 3 ), (13) 

where AH 2 = (T 0 |(H — £c) 2 |’4 , o) with E 0 = (’Fo|Hl’I , o)> 
the coefficient of t 2 in (12) can be interpreted as the 
variance of the energy in the initial state |d'o), A H = 

x /N[l + A(N-l)]/(2 V / 2). 

As a result of the free-expansion dynamics, the expo¬ 
nential regime [ ] is absent due to the lack of resonant 
states and a transition to the long-time behavior follows, 
when the survival probability is given by, 

/2 \ n[i+A(n-i)] 

5n,a(*) - (jJ ■ (14) 

Hence, the survival probability decays as a power-law in 
time. The power law exponent depends linearly on the 
interaction strength A (the GES parameter) and exhibits 
an at most quadratic dependence on the particle number 



OJ 0 t 


FIG. 1. Exact decay dynamics of a Calogero- 
Sutherland gas. The early stage exhibits a quadratic depen¬ 
dence on time and is followed by a smooth transition to the 
long-time power-law scaling. The dependence of the power- 
law exponent on the particle number (N = 2) changes from 
linear to quadratic as the interaction strength A is tuned 
from the non-interacting case to the Tonks-Girardeau limit 
(A = 1) and other values A ^ {0,1}, describing Haldane 
anyons with fractional exclusion statistics. From top to bot¬ 
tom, A = 0,1/2,1,3/2, 2. 


N. This is the main result of this manuscript. Its deriva¬ 
tion required (i) the scaling dynamics of the exact time- 
dependent coherent states, (ii) the use of the Bijl-Jastrow 
form of the initial state, and (iii) the identification of the 
leading term at long times. The scaling dynamics holds 
exactly for the CS gas and simplifies the ensuing analy¬ 
sis by contrast to other many-body systems such as, e.g., 
the ID Bose gas, where only recently moderate progress 
accounting for its dynamics has been reported [52-55]. 
Comparing the first leading terms in a long-time asymp¬ 
totic expansion, it is found that the power-law (14) sets 
in when the time of evolution satisfies 

t. > v / 2N(1 + A(N-1)). (15) 

For free-bosons, the power-law exponent becomes lin¬ 
ear in the particle number 

S N ,o(t) ~ (2A) N - ( 16 ) 

The case of hard-core bosons correspond to A = 1, and 
leads to a power-law exponent quadratic in the particle 
number 

S n ,i(£)~( 2/£) n2 . (17) 

The change in the scaling with N was conjectured an¬ 
alyzing quasi-free systems of N = 2, 3 particles [16, 17, 
22, 23, 25]. Eqs. (14), (16), and (17) prove that this is 
indeed the case. 
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III. ROBUSTNESS OF THE SCALING 

It is worth emphasizing that the power-law behavior 
(14) is observed in the long-time dynamics of other multi¬ 
particle observables such as the non-escape probability 
from a region of space, e.g., where the initial state is 
initially localized. Explicitly, we define the N-particle 
non-escape probability as 

N 

Pn,aW := / TTI^^ > O'N; *) 1 2 , (18) 

„=i 

which is the probability for the N particles to be found si¬ 
multaneously in the A-region and can be extracted from 
the full-counting statistics [23]. Explicit computation 
shows that 

PN A(*) = t N[l+A(N-l)] J NA(*)- ( 19 ) 

Taking A = [—a/2,a/2], the integral /n,a(£) becomes 
time-independent at long expansion times when b 
\/Aa, i.e. t |Aa 2 — ll 1 / 2 , 

N „? 

\\&qie~^ \ Qi - qj \ 2X , 

i—1 j>i 

~ 0 n[i+a(n-i)] 5n ( 1)1iA ) ) ( 20 ) 

where S n (a, /3, 7 ) is the Selberg integral [49]. As a result, 
the same power-law scaling sets in, i.e., 

^n,a (0 ex iSn.aW t _N[1+A(N_1)1 - ( 21 ) 

However, the dependence of the power-law exponent on 

N and A is lost when studying the decay in terms of one- 
body observables such as the one-particle density profile 
n(x,t) integrated over the region of interest A, 

P{t)~ [ n(q,t)dq. (22) 

J A 

To illustrate this, let us consider the exact evolution 
of the density profile that under scaling dynamics is 
given by n(q,t) = n(q/b,0)/b. Although an explicit 
computation of the density profile n(q, 0 ) is possible in 
the CS model, it would suffice to consider the large N 
limit. Then, n(q, 0) follows Wigner’s semicircular dis¬ 
tribution n(q, 0 ) = \J 1 — g 2 2 N/ 7 r which is already inde¬ 
pendent of A. Under free expansion, it is found that 
pit) ~ 2 aN/( 7 rf), where the power-law exponent is inde¬ 
pendent of N. The same conclusion holds when using the 
expressions for low N and large A available in the litera¬ 
ture for n(q, 0) [49, 56], but for the fact that the prefactor 
acquires a dependence on A/N. Generally, the 1/t power- 
law decay of pit) can be expected as the density profile 
flattens out at long expansion times, becoming approxi¬ 
mately constant over the region A, so that the integrated 
density profile pit) is governed by the normalization fac¬ 
tor 1 /bit). 





ui 0 t w B r 



uj 0 t 
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FIG. 2. Multi-particle state reconstruction. The sur¬ 
vival probability at a given evolution time is the sum of the 
product of the survival probabilities Sn.a^ — t)<Sn,a(t) (or¬ 
ange), the memory term M{t,r) (green) and the interference 
term I{t,r) (red), as it follows from Eq. (24). All contribu¬ 
tions are normalized to normalized to the value at the final 
time iSn,a(I) and the blue line is set at unity to identify the 
dominant term, (a) During the decay of a single particle all 
terms are significant, (b) As the particle number is increased, 
the memory term becomes dominant (N = 3, A = 1). (c) For 
larger values of the interaction strength, state reconstruction 
describes accurately the decay dynamics except for values of 
r which are small or comparable to the total evolution time 
t (N = 3, A = 2). (d) Increasing the particle number (N = 6, 
A = 2) further prolongs the interval governed by state recon¬ 
struction. The dynamics is induced by suddenly switching off 
a harmonic trap with initial frequency cjq, t = 15/wq- 


One might also wonder whether a non-sudden modu¬ 
lation of the trapping frequency will affect the long time 
power-law behavior. We show next that as long as the 
frequency of the trap is permanently switched off after a 
given time t = to, the same power law scaling sets in. In¬ 
deed, assume that b(t 0 ) = b 0 , b(t 0 ) = v 0 , K[t > t 0 ) = 0 . 
Then, 


b(t) = 


(60 + v 0 (t — t 0 )) 2 + 


(t-t 0 y- 


b 2 0 


(23) 


which tends to bit) ~ t(vg + 1 /frg ) 1 / 2 at large expansion 
times. As a result, only the prefactors of the survival and 
nonescape probability are affected, and the scaling is still 
dictated by ( 21 ). 


IV. MANY-PARTICLE STATE 
RECONSTRUCTION 

We next analyze the relevance of state reconstruction 
in the CS model as an example of a many-particle sys¬ 
tem. Using the Ersak equation (3) [9, 33], the following 
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decomposition of survival probability is obtained 

Sn.aW = S n ,a H ~ +X(t,T), 


where the first two terms admit a classical interpreta¬ 
tion. In particular, <Sn,a(* — t)«Sn,a(i") is the proba¬ 
bility for the system to survive in the initial state at 
time t provided that it was in the initial state at time 
r. Similarly, = \M(t,T)\ 1 2 3 4 5 6 7 8 9 10 11 12 accounts for state 

reconstruction in a classical sense, i.e., it is the prob¬ 
ability that the state has decayed at time r and re¬ 
construct the initial state at time t. The last term in 
(24) represents the interference between the amplitudes 
for the two histories just described, i.e., In,a( t, t) = 
23?[M(i, — t)^4n,a(t)]- Here, the survival am¬ 
plitude is 


*4n,a (t) 



ib)- 


-ir(t) I f[l+A(N-l)] 


(24) 


Figure 2 analyzes the relevance of each term in (24) nor¬ 
malized to <Sn,a(£) and after gauging away the dynami¬ 
cal phase EoT{t)/h = IV [1 + A (IV — l)]r(t)/2. Different 
regimes can be distinguished as a function of the param¬ 
eter /? = N[1 + A(N — 1)]. For a fixed evolution time t and 
small values of /3 all terms in (24) play a role. For larger 
values of /?, achievable by increasing either the particle 
number or the interaction strength, A4(f, r) dominates 
the contribution to the survival probability. Thus, the 
state reconstruction governs the long-time decay, except 
for values of t/r close to {0,1}, when all processes remain 
relevant. 


In conclusion, we have characterized the the exact de¬ 
cay of an interacting many-body quantum fluid released 
from a harmonic trap. Exploiting the self-similarity of 
the ensuing dynamics, the long-time power-law behav¬ 
ior of the survival probability was shown to be dictated 
by the strength of the interactions even at arbitrarily 
large expansion times. The scaling of the power-law ex¬ 
ponent is at most quadratic in the particle number. The 
non-exponential character of the decay can be attributed 
to the many-particle state reconstruction of the initial 
state. Our results can be extended to other systems 
including SU(v) spin degrees of freedom and fermionic 
exchange statistics [39]. As an outlook, it is worth ex¬ 
ploring higher dimensional systems like the 2D Bose gas, 
for which self-similar dynamics holds [57] up to quantum 
anomalies [58] and the generalized exclusion parameter is 
known [59]. In systems lacking self-similar dynamics, the 
role of interactions can be disentangled from that of the 
exclusion statistics and further studies will be illuminat¬ 
ing. A prominent example is the one-dimensional Bose 
gas with contact interactions, where recent advances in 
describing its dynamics have been reported [52-55]. 
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